Triangular Spectral Element simulation of two-dimensional elastic wave propagation using unstructured triangular grids
نویسندگان
چکیده
S U M M A R Y A Triangular Spectral Element Method (TSEM) is presented to simulate elastic wave propagation using an unstructured triangulation of the physical domain. TSEM makes use of a variational formulation of elastodynamics based on unstructured straight-sided triangles that allow enhanced flexibility when dealing with complex geometries and velocity structures. The displacement field is expanded into a high-order polynomial spectral approximation over each triangular subdomain. Continuity between the subdomains of the triangulation is enforced using a multidimensional Lagrangian interpolation built on a set of Fekete points of the triangle. High-order accuracy is achieved by resorting to an analytical computation of the associated internal product and bilinear forms leading to a non-diagonal mass matrix formulation. Therefore, the time stepping involves the solution of a sparse linear algebraic system even in the explicit case. In this paper the accuracy and the geometrical flexibility of the TSEM is explored. Comparison with classical spectral elements on quadrangular grids shows similar results in terms of accuracy and stability even for long simulations. Surface and interface waves are shown to be accurately modelled even in the case of complex topography with the TSEM. Numerical results are presented for 2-D canonical examples as well as more specific problems, such as 2-D elastic wave scattering by a cylinder embedded in an homogeneous half-space. They all illustrate the enhanced geometrical flexibility of the TSEM. This clearly suggests the need of further investigations in computational seismology specifically targeted towards efficient implementations of the TSEM both in the time and the frequency domain.
منابع مشابه
AIAA 2001–2623 Multi-Dimensional Upwind Constrained Transport on Unstructured Grids for ‘Shallow Water’ Magnetohydrodynamics
Novel Multi-dimensional Upwind Constrained Transport (MUCT) schemes on un-structured triangular grids are described. Constrained Transport (CT) discretizations conserve the divergence-free nature of divergence-free vector fields on the discrete level. Multi-dimensional Upwind (MU) schemes generalize the concept of dimensionally split upwind schemes for hyperbolic systems to truly multidimension...
متن کاملStress Waves in a Generalized Thermo Elastic Polygonal Plate of Inner and Outer Cross Sections
The stress wave propagation in a generalized thermoelastic polygonal plate of inner and outer cross sections is studied using the Fourier expansion collocation method. The wave equation of motion based on two-dimensional theory of elasticity is applied under the plane strain assumption of generalized thermoelastic plate of polygonal shape, composed of homogeneous isotropic material. The freque...
متن کاملA nite element method for the two-dimensional extended Boussinesq equations
A new numerical method for Nwogu’s (ASCE Journal of Waterway, Port, Coastal and Ocean Engineering 1993; 119:618) two-dimensional extended Boussinesq equations is presented using a linear triangular nite element spatial discretization coupled with a sophisticated adaptive time integration package. The authors have previously presented a nite element method for the one-dimensional form of these e...
متن کاملInternational Journal for Numerical Methods in Fluids
The weak Lagrange–Galerkin finite element method for the two-dimensional shallow water equations on adaptive unstructured grids is presented. The equations are written in conservation form and the domains are discretized using triangular elements. Lagrangian methods integrate the governing equations along the characteristic curves, thus being well suited for resolving the non-linearities introd...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2006